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Phase estimates in adaptive-optics systems are computed by use of wavefront 
sensors such as Shack-Hartmann or curvature sensors. In either case the 
standard error of the phase estimates is proportional to the standard error 
of the measurements; but the error-propagation factors are different. We 
calculate the ratio of these factors for curvature and Shack-Hartmann sensors 
in dependence on the number of sensors, n, on a circular aperture. If the 
sensor spacing is kept constant and the pupil is enlarged, the ratio increases as 
n*^'^. When more sensing elements are accommodated on the same aperture, it 
increases even faster, viz. proportional to n*^'^. With large numbers of sensing 
elements this increase can limit the applicability of curvature sensors, (c) 
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1. Introduction 

In adaptive-optics (AO) systems the wavefront sensing is commonly performed with Shack- 
Hartmann (SH) sensors [1] or Roddier curvature sensors (RC) [2]. Both types of wavefront 
sensors have been successfully employed in fields such as astronomy [3] or ophthalmology [4] , 
[5], [6] - with up to 188 sensing elements in the case of RC sensors [7] and nearly a thousand 
elements in the case of SH sensors [8] . For the next generation of adaptive optic systems, with 
increasingly larger numbers of elements, it is of interest to compare - for different numbers, 
n, of sensing elements - the two techniques in terms of the error-propagation factor, i.e. the 
ratio of the standard error of the inferred phase estimates to that of the measured signal. 
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The error-propagation factor, csh; for SH sensors has been found to increase logarithmically 
with the number, n, of sensing elements when their spacing is kept constant and the pupil 
surface is increased [9], [10]. Simulations by N. Roddier have shown that under the same 
condition, the error propagation, crc; for R-C sensors increases linearly with n [11]. For a 
constant pupil surface - i.e. if more sensing elements are packed on the same aperture - crc 
is almost independent of n [12]. 

Here, the ratio, e^'-^ /e^^ , of the error-propagation factors for RC and SH sensors is exam- 
ined in relation to the number of sensing elements, on a circular aperture. The ratio is found 
to increase as nP'^ when the spacing of the sensing elements is fixed. It increases even faster, 
roughly as n^'^, when more sensing elements are placed on a given aperture. 

2. General considerations 

In practice the wavefront reconstruction is performed in terms of matrices that are pre- 
calculated for the particular AO system. A command matrix of voltages is sent to the set of 
mirror actuators so that in combination with their influence function, i.e. the matrix of phase 
shifts caused by a unit shift (voltage increment) of an actuator, the correct reconstruction 
matrix of phase increments is obtained. In other words, the reconstruction matrix that 
transforms the measured slope or curvature values into the phase estimates is realized as the 
'product' of the command matrix and the influence function. 

While the actual reconstruction algorithms differ, they all contain - regardless whether 
they are used with SH- or RC-sensors - the solution of the Laplace equation in terms of 
mean- curvature values. Separation of the error propagation into the two inherent steps sensor 
signal to curvature and then curvature to phase estimates will, thus, make the analysis more 
transparent without changing the results. 

For an orthogonal grid of sensors the linkage between slope and curvature values is straight- 
forward. For the common radially symmetric arrangements of curvature sensors it is more 
complicated. The reconstruction formalism is then still based on the solution of the Laplace 
equation in terms of the mean curvatures, and the algorithm that is described in section 2.C 
(and is employed in sections 4 and 5) will still apply in just slightly modified form. However 
the separation of the two steps of error propagation is most transparent for an orthogonal 
pattern of sensors, and this special case will bring out the essential characteristics of error 
propagation for SH and for RC sensors. 

2. A. Phases, slopes, curvatures: 

Let the phases at the sensor position {i, k) be Fj fc. Only the phase differences are of interest. 
They are measured in x- and y-direction by SH sensors and are termed slopes. The relation 
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between the true, i.e. error-free values are: 

•^i,k Fi+l,k Fi,k 

yi,k = Fi,k+1 — Fi,k (1) 

The mean curvatures are measured by RC sensors: 

Ci,k = {xi^k - Xi^i^k + yi,k - yi,fc-i)/4 = + + Fi^k~i + -Fi,A.-+i)/4 - -Fi./t (2) 

Analogous relations that include less than 4 adjacent sensors apply to the sensors at the 
boundary of the aperture. 

In actuality, the variables Xi^k or i/i^k are slopes times grid spacing, d, and the are mean 
curvatures times (P, i.e. all variables are phase increments. To distinguish the variables slope 
and curvature from the slope and curvature proper they are written in italics, the latter are 
written in roman. 

2.B. Iterative reconstruction 

As shown in earlier analyses (e.g. Hudgin [9]), each phase estimate, Fi^k, is a linear combi- 
nation of the elements xi^m and yi^m of the slope matrices: 

Fi,k = ^ ^/i,fc(^ m) xi^rn + ^ ^fi,k{^, m) yi^rn (3) 

l,m l,m 

^ fi^k and ^ fi^k are the reconstruction matrices from slopes to phase estimates, fc. With 
regard to the mean curvatures one can use the same concept, but with different reconstruction 
matrices : 

Fi,k = ^ fi,kil, m) Cl^rn (4) 
l,m 

Even where slopes are measured they are converted into curvatures , because the less com- 
plicated Eq. 4 provides the same phase estimates as Eq. 3. The reconstruction matrices f\k 
depend on the sensor grid and position; they are pre-calculated for a given AO setup. 

A given set of curvatures determines the best phase estimates, i.e. the estimates that are 
the least-squares fit to the measured curvatures or the underlying slopes. Starting with an 
arbitrary set, Fj fc, of estimates the value on every grid point is replaced in each iteration 
step by the average of the 4 adjacent values, with adjustment for the measured slope or the 
corresponding curvature data [9]: 

Fi^k ^ ( Fi^i^k + Fi^i^k + Fi k-i + Fi k+i + Xi,k — Xi-i^k + yi,k — l/i,fc-i)/4 

= (Fi_i,fc + Fi^i^k + Fi^k-i + Fi^k+i)/4: — Ci^k (5) 
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The new value Fj ^ fits best the preceding phase estimates at the adjacent 4 grid points 
and the measured slope or curvature data. The analogous relations for sensor positions on 
the edge of the pupil which have only 1, 2, or 3 neighbors are evident and are not listed. The 
interior curvature values contain no information on tip and tilt. The procedure converges to 
the least-squares solution. Depending on n, several hundred iterations may be required, but 
the procedure converges correctly. 

2. C. Reconstruction by matrix inversion 

The iterative procedure is simple and transparent, but the direct solution by inversion of 
one large matrix that connects phases and mean curvatures is a convenient alternative. To 
employ it one arranges the variables for the different sensor positions as vectors, rather than 
matrices. This is merely a change of notation; instead of dealing with a matrix of matrices 
as in Eq. 4 one deals with a single matrix equation: 



In this alternative description F and c are the vectors of the phase values and the mean 
curvatures. They contain n elements, where n is the number of sensors, x and y is the 
analogous notation for the slopes. The nxn matrix U contains in each row at the appropriate 
positions the term -1 and four terms 1/4 when the phase corresponds to an interior position. 
For positions on the boundary of the sensor grid there are, apart from -1, three terms 1/3, 
two terms 1/2, or one term 1. 

If all phases are equal, the curvature vector, c, is zero, which shows that the matrix U is 
singular. The reconstruction relation contains, therefore, instead of the familiar inverse the 
pseudo-inverse, R, of U: 



Eq. 7 is the equivalent of Eq. 4. The column of the matrix R that corresponds to the sensor 
position (/, m) in Eq. 4 is a vectorial representation of the matrix m). Also it is readily 
seen that the equivalent of Eq. 3 is : 



where each column in the two matrices ^R and ^R is the difference of the two columns in 
R that correspond to the pair of sensor positions associated with the slope x or y. 

The vectorial representation of the slopes, mean curvatures and phases is utilized only in 
section 5, where Eqs. 7 and 8 are used to obtain, for sensor numbers up to n = 500, the 



c=U-F 



(6) 



F = R-c 



(7) 



F = ""R-x+yR-y 
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error-propagation factors for RC and SH sensors. The use there of the inversion algorithm is 
a matter of convenience. In the absence of tools to handle very large matrices the iterative 
method serves the same purpose. 

2.D. Error propagation 

The different types of wavefront sensors have particular strengths and limitations that de- 
termine their applicability. The present discussion is directed at one characteristic of SH and 
RC sensors: the different propagation of measurement errors to phase-estimation errors. 

Let z and v be the measured signals of the SH and RC sensors; z and v are dimensionless 
values that can - with coefficients specified in the subsequent section - be converted to the 
slope and mean curvature values. 

If az and are the standard errors of the SH and RC signals, z and v, then the phase 
errors, (Jf,sh and (Jf,rci are the products of the errors and and the error-propagation 
factors, csH and crc- 

0'F,SH = O'z ■ esH 

(^F,RC = 0"d ■ GRC (9) 

The exploration of these error-propagation factors will, as indicated in Fig. 1, be separated 
into two steps. The first step considers the two fairly trivial partial factors from and cx^ 
to the standard errors of the mean curvatures: 

ei,Rc = crc,Rc/<yv (10) 

This step requires the specification of the calibration factors from SH- or RC-sensor signal 
to mean curvature, which will be detailed in the next section. 

The second step - which is the major part of this study - is the determination of the 
propagation factors from the mean curvature errors to the phase errors, crp^sH and crpjic'. 

^2,SH = o'f,sh/o'c,sh 

e2,RC = (^F,RC / (^c,RC (H) 

The error propagation for RC and for SH sensors will be compared in terms of the ratio: 

r = enc/esH = ei,/?c'/ei,sif ■ e2,Rc/e2,sH = n ■ r2 (12) 

^i,RC, ^i,SH and ri are considered in the next section in relation to the number, n, of 
sensors. Sections 4 and 5 will subsequently deal with the second step, the error propagation 
from mean curvature to phase. 
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3. Comparison of the curvature errors from SH and RC sensors 

3. A. SH sensors 

The signal, z, equals the measured displacement of the image in the detector plane - in units 
of pixel size, p. The relation between the slope value s, i.e. Xi^k or ?/j fc, and the signal, z, of 
the SH sensor is [1]: 

. = (13) 

where A is the wavelength, fi the focal distance of the SH lenslets, and d the grid spacing. 

The same relation applies to the standard errors, cr^ and a^, of s and z. To obtain the 
expression for ci^sh one needs, in addition, the relation between as and the standard error, 
o'cSHi of the mean curvature derived from the slope values. This relation is simple, since the 
standard error of the mean curvature equals - for internal sensors - half the standard error 
of the slope values (see Eq. 2). For the sensors on the boundary the value is somewhat larger, 
but averaged over all sensors the difference can be disregarded - especially for large n. With 
the factor 1/2 one obtains: 



ei,sH = (TcSh/cts = ndp/Xfi (14) 

when n is increased, while the grid distance, d, and the size of the sub-apertures remain the 
same. A, p and // are kept constant. The error-propagation factor, ci^sh, is then independent 
of n. 

On the other hand, when more sensors are packed on the same aperture, the focal length, 
//, of the Shack-Hartmann lenses is generally increased in order to keep the amplitude of 
the image motions constant in the detector plane. Since the amplitude of the sensed phase 
distortions decreases with sub-aperture size as ~ n ^'^ [13], // is adjusted proportional to 
j^5/i2^ One concludes that the error propagation from the SH sensor signal to mean curvature 
is then roughly inverse to the number of sensors: 

ei,sH oc n-''/'' (15) 

3.B. Curvature sensors 

The relation between the mean curvature, c, and the signal, v, of the RC sensors is [2]: 

27r d^l 

Accordingly: 

27r d^l 
T 7^ 



c=——v (16) 



ei,Rc = crc,Rc/o-v = — (17) 
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Again, the error-propagation factor is independent of n, when n is increased, while the 
grid distance, d, and the size of the sub-apertures remain the same. 

When more sensors are packed on the same aperture, the telescope focal length, /, remains 
unchanged, but the extra-focal distance, /, of the curvature sensor is increased proportional 
to n^/^, so that the size dl/ f, of the sensing elements in the detector plane stays constant. 
It follows, that 

Eqs. 15 and 18 imply that the ratio, ri, increases when more sensors, n, are packed on the 
same aperture: 

ri = ei,ijc/ei,5if ~ n^^^^ (19) 

4. Error propagation from curvature to phase — Comparison of SH and RC 
sensors 

It is not surprising that the error propagation from curvature to phase estimation should be 
somewhat worse for the directly determined mean curvatures . Broadly speaking, a measure- 
ment error of the slope at a sensor position causes an imprecise estimate of the local phase 
increment only, while an error of the measured curvature causes an imprecise estimate of the 
rate of phase increments, which then affects also distant estimates of phase increments. This 
is analogous to reconstructing a function of a single variable from its imprecisely measured 
second derivative, which is less stable than reconstructing it from the imprecisely meas- 
ured first derivative. While there is, thus, no particular need to explain the difference, it is 
still helpful to consider the issue in analytical terms before numerical results are derived in 
sections. 

4-. A. Correlated vs. uncorrelated curvatures 

The observed slopes or mean curvatures and the phase estimates are the sum of the true 
values and the errors. Since the reconstruction is additive, the errors of the phase estimates 
depend - regardless of the true slope, curvature and phase values - only on the measurement 
errors. The subsequent considerations are therefore simplified by taking the true values to 
be zero, so that the slope and mean curvature values Xi^mi yi,m and stand for the errors 
only. 

To facilitate the comparison of error propagation, two different procedures will be invoked 
that generate sets of mean curvatures as obtained with SH and RC sensors. Both procedures 
start out by simulating independent and unbiased random slope errors; the mean square of 
these errors is set equal to 4, so that the mean square of the corresponding ci^m is unity. The 
boundary effect is here disregarded, i.e. all measured curvatures are treated as interior values 
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that are linked to 4 slope values. The boundary could be accounted for, but this aspect is 
here not essential. 

Procedure A corresponds to the use of SH sensors; it generates one set of 2n slope errors, 
Xi,m, yi,mi and derives from these values, in line with Eq. 2, the set of mean curvatures , ci^m- 
Adjacent values ci^m are then to some degree correlated. 

Procedure B simulates separately for each position, {l,m), the four associated random 
errors Xi_i^rn,yi-i,m,xi^rn, and yi^m, and derives from them the curvature error, • As with 
procedure A, the resulting mean square curvature error, < >, is unity, but the are 
now uncorrelated, i.e. they resemble a set of mean curvatures obtained with RC sensors. 

Consider first the case of merely one single slope error, say Xn,m = With procedure A this 
error, Xn,m = e, corresponds to the two mean curvature errors, Cn,m = e/4 and c„+i_m = — e/4. 
The addition of these two terms according to Eq. 4 gives the associated phase errors: 

= e/4 ■ (/,,,(/, m) - /,,,(/ + 1, m)) (20) 

One obtains the corresponding mean square phase error: 

dl^ = eVl6-[<(/(/,m)-/(/ + l,m))2> -<(/(/, m)-/(/ + l,m))>2] 

= eVl6 ■ [< /(/, m)2 > - < /(/, m) >2 + < /(/ + 1, mf >-</(/ + 1, m) >2 
-2(< /(/, m) . /(/ + 1, m) > - < /(/, m) >•</(/ + 1, m) >)] (21) 

< .. > denotes the average over all {i, k). 

With the notation var{) and cov{, ) for variance and covariance and with < >= 1, 
this takes the form: 

dh = [var{f{l, m)) + var{f{l + 1, m)) - 2cov{f{l, m), /(/ + 1, m))]/4 (22) 

The total mean square phase error, <Tpg^, is obtained by summing Eq. 22 for all 2n slopes. 
The contributions of the slopes xi^m and yi^m being equal, it is sufficient to sum Eq. 22 over 
all n positions (/, m) and then multiply the result by 2 : 

^F,SH = 4,SH = Yl^^^^^f^^' ~ cov{f{l, m), /(/ + 1, m)] (23) 

l,m 

o'F,SH equals the squared error-propagation factor e2,sH, because ac,sH = 1- 
With procedure B the single slope error Xn,m corresponds to two curvature errors, = 
ei/4 and Q+i^m = —^2/4. Both for ei and €2 the mean square error is < >, but the mixed 
term < ei ■ ei > in the equation below vanishes because of the statistical independence of ei 
and 62- The increment of the mean square phase error is thus: 

die = ■ [< fil, mf>-< fil, m) >''] + e?/16 ■[</(/ + 1, m)^ > - < /(/ + 1, m) 

-ei 62/8 ■ [< /(/, m) ■ /(/ + 1, m) > - < /(/, m) >■</(/ + 1, m) >] (24) 
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Since the mixed terms disappear one obtains instead of Eq. 23: 

<^F,RC = 4,RC = Yl ^'«'"(/(^' "^)) (25) 

l,m 

This is a notable result: The error-propagation factors differ by the covariance term in 
Eq. 23. This term converges with increasing n towards the variance, which keeps down the 
error propagation, when the number of SH sensing elements becomes larger. In the case of 
independent curvatures the covariance term is absent. This increases the error propagation 
with RC detectors and makes it roughly proportional to n, i.e. to the number of terms in 
Eq. 25. 

4.B. Analytical approximation 

The above conclusion, the poorer error propagation for curvature measurements with large 
sensor numbers, will now be quantified. For this purpose a large grid of sensors is considered 
that extends beyond a central region of n sensor positions for which the phase estimates 
are actually sought. This large-grid approximation is informative because it entails a simple 
reconstruction function, fi^k, that - regardless of the sensor position, (/,m) - is roughly 
proportional to the logarithm of the distance between sensor positions. 

While the large-grid approximation brings out the characteristic difference between the 
use of curvatures from SH and RC detectors, it provides results that differ somewhat from 
the more accurate numerical evaluation in section 5. 

With grid spacing d the distance of position (i, k) from the central position, (0, 0), is: 

dVfTl^ = d-u (26) 

The function f{x, y) = a ln(n) + b has mean curvature zero for u > 0, i.e. its Laplacian - 
or twice the mean curvature - vanishes: Af = - ^(u ^) = 0. 

u au ^ au ' 

f{x,y) is, therefore - at least at somewhat larger u, where the difference quotient equals 
the differential quotient - a solution for the reconstruction function. The constant term, b, is 
irrelevant and can be disregarded. This is in line with the observation that the reconstruction 
function is independent of d, which means that for given curvature errors the phase errors 
are the same, regardless whether the grid spacing, d, is decreased in order to accommodate 
a larger number of sensors, or whether d is kept constant and the aperture is enlarged. 

Numerical evaluation in terms of the iterative method shows that: 

= 2/7r ln(n) = \n{z^ + e)/TT and /o,o = /(O) = -1 (27) 

applies exactly at u = 0, u = 1 and for large u. For intermediate u{> 1) the difference to the 
exact solution is too small to be of concern for the present approximation. 



9 



Eq. 27 gives the reconstruction function, /j fc(0,0), for the central sensor position. For 
position (/,m) it is shifted accordingly: 

/,,fc(Z,m) = ln((z-/)2 + (A;-m)2)/7r and m) = -1 (28) 

Fig. 2 gives the error-propagation factors e2,Rc and 62,5// that are obtained by inserting 
the above expression into Eqs. 23 and 25, where the variance and covariance are - again with 
the symbol < > for average over all {i, k): 

var{f{l,m)) = < fi,k{l,rnf > - < fi,k{^,rn) 
cov{f{l,m)J{l + l,m)) = < /i,fc(/,m) ■ /i,fc(/ + l,m) > 

- < /,,fc(/,m) > • < /i,fc(/ + l,m) > (29) 

The results are given for symmetric circular sets of sensors on an orthogonal grid. The 
center of the grid is on the midpoint between 4 sensors ( n = 12,16,24,32,..). The main 
finding is the very slow increase with n of the error propagation from curvature to phase in 
the case of the indirectly measured curvatures (SH), versus the more substantial increase for 
the directly measured curvatures (RC). 

Fig. 3 gives the corresponding ratio r2 = e2,_Rc / G2,sh of the error-propagation factors. For 
sensor numbers up to about 40, r2 is less than 2, for n = 500 it is roughly 5.5. 

5. Explicit computation of the error propagation from curvature to phase 

Eqs. 7 and 8 can, largely in analogy to the preceding section, be used to derive exact values 
of the error-propagation factors e2,ijc and 62, 5^- Since the mean curvature errors from RC 
detectors are statistically independent, their mean-square contributions to the phase errors 
add up. The contribution of the sensor labeled k in the vector representation equals therefore 
the product of the mean-square curvature error < > and the sum of the squared com- 
ponents of column k in the reconstruction matrix, R. Since the contributions of all sensors 
add up, the mean squared phase error equals the sum of all squared components of R. The 
error-propagation factor is thus: 



e2,RC = n'-' < Rlu > (30) 

To compute the error propagation with SH detectors one must go back - as in the con- 
siderations in the previous section - to the slope errors, i.e. one needs to use Eq. 8 instead 
of Eq. 7. Each of the columns of the matrices ^R and ^R are, as pointed out in section 2. C, 
the differences of two associated columns of R. The sum of the squared matrix elements is, 
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therefore, in analogy to Eq. 23, the sum of the squares of the elements of R minus a sum of 
the product of associated elements. This provides the relation: 

e2,SH = n'-' (< ""Rl > + < ^4,. >)/4 

= n'-' < -Rl > /2 (31) 

The diamonds in Fig. 4 give the results. While the dependences on n agree essentially with 
those in Fig. 2, the absolute values of e2,Rc and 62,5'^ are larger than in Fig. 2. This is not 
surprising, because the approximation in section 4. B has inherent symmetry that excludes 
tip and tilt. For a further comparison it is therefore informative to determine the phase errors 
with tip and tilt correction. They are given in Fig. 4 by the smaller dots. The correction for 
tip and tilt reduces the phase-estimation errors substantially, especially for the RC data. 

The contribution of tip and tilt is seen to almost double the standard error of the phase 
estimates in the case of directly measured curvatures, which means that it contributes pre- 
dominantly to the variance of the phase-estimation errors. For curvature values derived in 
terms of SH sensors the relative contribution of tip and tilt is markedly less; it contributes 
about half of the variance. 

The analytical approximation in section 4. B has been solely based on the mean curvatures . 
Since these carry no information on tip and tilt, the curves in Fig. 2 should correspond to 
the results in Fig. 4 without tip and tilt. They run, in fact, below the curves that include tip 
and tilt, but lie above those without tip and tilt, which reflects inaccuracies due to the large 
grid- approximat ion . 

Overall, however, there is reasonable agreement, and the essential finding is, that the 
curves for the RC detectors in Figs. 2 and 4 have the same slope of about 0.5, which means 
that e2,Rc is proportional to n^'^. The curves for the SH detectors exhibit a slope close to 
0.1, i.e. they show that e2,Rc increases roughly proportional to over the relevant range 
of n the slow increase is equally in accord with a logarithmic dependence as earlier reported 
by Hudgin [9] and by Fried [10]. 

Fig. 5 gives the ratios r2 that corresponds to Fig. 4. With tip and tilt excluded r2 reaches 
the value 2 at = 50 and exceeds 4 at n = 500. The values without tip and tilt correction 
are about 1.5 times larger, but they are not highly relevant, because in AO the use of RC 
detectors is commonly combined with measurements to allow tip and tilt correction. 

6. Conclusion 

Shack-Hartmann slope detectors and Roddier curvature sensors have particular strengths 
and limitations that determine their applicability. The present analysis has been concerned 
with one characteristic of the wavefront sensors, the propagation of measurement errors to 
phase-estimation errors. With either type of sensor the standard error of the phase estimates 
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is proportional to the standard error of the curvature values, but the error-propagation fac- 
tors differ. The errors of the directly measured curvatures are statistically independent, and 
therefore their contributions to the mean-square phase errors add up. In contrast, when 
the mean curvature values are obtained from slope measurements, adjacent values are neg- 
atively correlated and accordingly their mean-square phase-estimation errors combine less 
than additively. 

The results confirm and quantify the earlier findings that the error propagation is larger 
when the curvatures are measured directly instead of being derived from measured slopes. 
The ratio, crc/^sh, of the error-propagation factors for curvature and slope measurements 
is seen to increase roughly as ^/n, when the size of the sensing elements is fixed and the 
pupil size is enlarged. When more sensing elements are packed onto the same pupil, the ratio 
increases even faster; Crc/^sh is then roughly proportional to n. 

Finally it needs to be noted that the standard errors of the phase estimates are the 
product of the standard errors of the detector response and the error-propagation factor (see 
Eq. 9). For photon-noise limited measurements, the standard error of the detector response 
is inversely proportional to the square root of the mean number of detected photons, i.e. the 
product of fluence and the area of the sub-apertures. Since the photon number is inversely 
proportional to n, one concludes that the standard errors of the phase estimates increase 
even faster with n than the error-propagation factors. 

The difference in error propagation is not decisive at small or moderate sensor numbers, 
n, and, generally, the error propagation is by no means the only or even the dominant aspect 
of the applicability of the different types of sensors. The substantial increase of the error- 
propagation factor for curvature detectors at large n can nevertheless be a critical factor 
that needs to be weighed against their various advantages. 
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Fig. 1. The analysis of the error propagation in SH- and RC-sensors is sepa- 
rated into two steps: signals to mean curvatures and mean curvatures to phase 
estimates. 
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Fig. 2. Upper line of symbols: e2,_Rc- Lower line of symbols: e2,sH- The results 
are given for the number n of sensors that occur within circular regions, with 
the grid-center positioned in between 4 sensors. 
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Fig. 3. Ratio r2 = e2,Rc / ^2,sh- The results are given for the number n of 
sensors that occur within circular regions, with the grid-center positioned in 
between 4 sensors. 
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Fig. 4. Error propagation factors 62, rc (upper pair of curves) and 62,5// (lower 
pair of curves). In each pair of curves the upper curve gives the total e2-factor, 
while the lower curve gives the reduced value that remains when tip and tilt 
are removed from the estimates. The results without tip and tilt compare to 
the ones in Fig.2. - The results are obtained for an orthogonal grid of n sensors 
within a circular aperture. The grid center is positioned in between 4 sensors. 
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Fig. 5. The ratio, r2, of the error-propagation factors e2,_Rc and e2^sH repre- 
sented in Fig. 4. The upper hne of symbols relates to the entire standard error 
of the phase estimates. The lower line of symbols relates to the standard error 
without the contribution of tip and tilt. It is in reasonable agreement with the 
approximate results represented in Fig. 3. 
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